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O ■ 

■ The Density Matrix Renormalization Group (DMRG) method with periodic boundary conditions 

is introduced for two dimensional classical spin models. It is shown that this method is more 
suitable for derivation of the properties of infinite 2D systems than the DMRG with open boundary 
conditions despite the latter describes much better strips of finite width. For calculation at criticality, 
phenomenological renormalization at finite strips is used together with a criterion for optimum strip 
\ width for a given order of approximation. For this width the critical temperature of 2D Ising model 

is estimated with seven-digit accuracy for not too large order of approximation. Similar precision 
(— I I is reached for critical indices. These results exceed the accuracy of similar calculations for DMRG 

. with open boundary conditions by several orders of magnitude. 
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^ . I. INTRODUCTION 

In 1992 the density matrix renormalization group (DMRG) technique in real space was invented by S. R. White 
101 and has been mostly applied to the diagonalization of one-dimensional (ID) quantum chain spin Hamiltonians. 
Three years later the DMRG method was redesigned by T. Nishino Q and applied to classical spin 2D models. The 
Ch DMRG method for the classical models is based on renormalization of the transfer matrix. It is a variational method 
O maximizing the partition function using a limited number of degrees of freedom and its variational state is written as 
I I ■ a product of local matrices . 

\ DMRG has been used for many various quantum models. It provides results with remarkable accuracy for larger 
^""^ ■ systems than it is possible to study using standard diagonalization methods. The 2D classical systems treated by the 
>^ ' DMRG method exceeds the classical Monte Carlo approach in accuracy, speed, and size of the systems j|] . A further 
. DMRG improvement of the classical systems is based on Baxter's corner transfer matrix the CTMRG 0, and its 
' generalization to any dimension jsj . 

, Applications of the DMRG technique for calculation of the thermodynamic properties in the 2D classical systems 
has been done by P,p|-pT[ . Treating of non-symmetric transfer matrices or non-hermitian quantum Hamiltonians has 
O '. also been studied by the DMRG technique 

• It was shown that DRMG method yields very accurate estimations of ground state energy of finite quantum chains 
and free energy of classical strips of finite width with open boundary conditions. We have developed the DMRG 
method with periodic boundary conditions for strips of classical spins and shown that, similarly as for quantum 
chains, it gives these quantities with much less degree of accuracy. Nevertheless, the DMRG method is mostly used 
for prediction of physical quantities and critical properties of infinite systems in connection with finite size scaling 
or extrapolation of the results from finite-size systems to infinite ones. The objective of this paper is to study the 
Q , DMRG and exact methods with two different boundary conditions for finite strips of various widths and compare 
O ■ their results with known exact results for infinite 2D system. It is shown that while for the exact diagonalization of 
^ I finite-strip transfer matrices scaling properties of the system improve, for DMRG approach there exists an optimum 
width for each degree of approximation. The developed approach is tested on 2D Ising model. 

Paper is organized as follows: in Sec. II we mention briefiy the DMRG for the open boundary conditions; Sec. Ill 
contains the modification of the DMRG method for the periodic boundary conditions; in Sec. IV we present the results 
obtained by the DMRG method with periodic and open boundary conditions, the exact diagonalization method and 
show how to determine the optimum strip width for finite-size scaling, and in Sec. V the results will be summarized. 



II. DMRG WITH OPEN BOUNDARY CONDITIONS 



The transfer matrix approach is a powerful method for exact numerical calculation of thcrmodynamical properties 
of lattice spin models defined on finite- width strips. If the width of the strip is too large and the capacity of the 
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computer is exceeded, the DMRG method is found to be useful for an effective reduction of the transfer matrix size. 
It can be used for calculation of global quantities such as free energy as well as of a spatial dependence across the 
strip of local quantities, e.g. spin correlation functions. 

The properties of an infinite strip of finite width L are given by the solution of 'left' eigenvectors and corresponding 
eigenvalues of the transfer matrix equation 

Y,^i{{a})T{{a}\W})^\^i{W}), (1) 

where {cr} is a set of L spins {cri,(T2, . . . defined on a row and {cr'} is a set of L spins on the adjacent row. 
The transfer matrix is a product of Boltzmann weights given by the lattice Hamiltonian. For non-symmetric transfer 
matrices besides the left eigenvectors the right eigenvectors 5'* should be calculated, as well. 
Reducing the size of the transfer matrix the standard DMRG technique proceeds in two regimes: 

(1) In the process of iterations, the infinite system method (ISM) pushes both ends of the transfer matrix further so 
that each step of the ISM enlarges the lattice size by two sites. The transfer matrix (superblock) is constructed from 
three blocks: left T; and right transfer matrices (blocks) and a Boltzmann weight Wb, in particular 

T[2,+2] =r/^VBTW, (2) 

where the index on the left hand side denotes the number of sites in one row of the whole superblock T at the jth 
step of iteration. The Boltzmann weight usually is a function of several spins interacting among each other, e.g. for 
the Ising model with nearest-neighbor interactions the Boltzmann weight has the form 

WB{cTia2\a[a'2) = "^^P |^'^~^('^i^2 + 0-^0-2 + (Tia[ + cr20-2)| ■ (3) 

In the first step of the ISM (for details, see [|l|,D) t/^^ = Tr^^ = Wb is put. The whole procedure has L/2 ~ 1 steps 

r/^^sTi'' ^ r/^VsTi^) ^ . . . -> T/^/^-iVsri^/"-'). (4) 

and stops when the desired strip width of L sites is reached. 

The first steps of the iteration scheme (^ are exact but if the superblock matrix T becomes too large, a reduction 
procedure, to keep the size of superblock constant, should be introduced. 

The first step of (4) introduces open conditions at the strip boundaries. If the temperature of the system is lower 
than the critical one and the strip width is wide enough, the symmetry of the system is spontaneously broken (order 
parameter becomes non-zero), and after reaching the fixed point of the iteration procedure, the system does not 
depend on the boundary conditions any more. The calculations with periodic boundary conditions described in the 
next Section give in this regime the same result as with the free ones. 

(2) the finite system method (FSM) improves numerical accuracy of ISM result by left and right moves (sweeps) 
according to the following prescription: 



t/^-'VbT.^i) ^ r/^-'VBr,(2) ^ . . . ^ r/^/'-^VsTi^/'-'^ (6) 

In the right sweep (5) the left blocks T; are calculated in the previous step of the sweep and the right blocks Tr are 
taken from the previous left sweep (in the first right sweep from ISM); similarly for the left sweep. 

The values of local thermodynamical quantities given by particular superblocks in the final sweep (after the 
steady state is reached) are spatially dependent. The values given by the superblock in the middle of the strip 

are the closest to the bulk ones. In this sense, the best transfer matrix eigenvalues as well as 
eigenvectors are those of the above-mentioned central superblock. The two largest eigenvalues are used for further 
finite-size scaling or extrapolation treatment. 



III. DMRG WITH PERIODIC BOUNDARY CONDITIONS 



The translational invariance of the infinite lattice is preserved in finite strips with periodic boundary conditions 
when strip boundaries are connected with bulk intersite interactions. In this case the strip forms an infinitely long 
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cylinder. If the radius of the cyhnder is small enough, the model can be easily solved by exact numerical diagonalization 
methods. 

In DMRG language, imposing periodic boundary conditions means that we have to introduce properly the connec- 
tion of both ends of the superblock transfer matrix T. Thus, in distinction to open-boundary case the superblock is 
constructed from two Boltzmann weights connecting two blocks at both ends (see Fig. 1, the rightmost diagram). 

xM/B(a,+4a,+3|a;.+4a^.+3)Ti^)(a,+3era2|a;-+3ea^)M^B(a2fTi|fT^al), (7) 

where the block spin variable — {1, 2, . . . ,to}, and the primed variables are denoted by filled circles and ovals in 
Fig. 1. 

In the first few steps the lattice is enlarged to the desired size; no degrees-of-freedom reduction is performed and 
the superblock transfer matrix remains equivalent to the exact one. As depicted in Fig. 1, the ISM starts with 



y(6) ^ Tl >WBTy-'WB defined on twelve sites where T, 
sites are added in each of further steps. 



(1) 



T, 



(1) 



WbWb, and one Boltzmann weight, i.e. four new 





FIG. 1. The first j steps of the ISM for the strip with the periodic boundary conditions. 



If 2^ > 771, the number of degrees of freedom should be reduced at each jth step to keep the order of the superblock 
matrix constant and equal to 2^* • . 

Summation in the equation for eigenvectors (1) of the transfer matrix (8) can be performed in two steps 

$(criCiCTj+4CTj+3Cr^+3C0'20'l) = ^ri^T j+zS.r(T2Wj+:i£.W2)W b{<J2(Ji\(jWi) 

X^'(cri6o-j+4CTi+3CrO-2) 

*K^,'a^.+4a;-+3C^^) = r'(^i^''^J+4ki^;^^+4)W^i3('^j+4fT,+3|a;-+4a;.+3) (8) 



which is depicted graphically in Fig. ^ 



^^3 



^j+3 



X $ (CTI CTj +4 CTj+3 cr' +3 ^^2 ) 




^j+3 




^4 



(5i ai 




^j+4 «Jj+3 



^jV4 ^V3 



FIG. 2. Graphical representation of Eq. The variables represented by filled circles and rectangles are summed over. The 
spins CTi and a'l at both ends of the superblock must be identified due to the periodic boundary conditions. 
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This procedure uses the left and right transfer matrix blocks to calculate properly the left and right eigenvectors 
and '^r, respectively, of the whole superblock for the periodic boundary conditions. Once we have the '^i and ^f^, 
the left and right density matrices can be constructed. 

Pl{£,lCrj+4\ClCrj+d = *i(c^l6o-J+40-j+3CrO-2)^'r(criCiCTj+4CTj+3^rCr2) (9) 

Pr(CTj+3^rkj+3C) = *i ('^l6o-J+4Crj+3CrO-2 )*r (o-l6CTj+4Crj+302 ) , (10) 

and by its complete diagonalization 

sets of left and right eigenvectors stored in Oi and Qi matrices, respectively, is obtained (analogously, for and 
Qr)- The indices i,j = 1, 2, . . . , 2m) run over all states of m-state multi-spin variable ^ and two-state variable 
a. For the last steps of ISM and all FSM steps half of the eigenvectors (corresponding to their lowest eigenvalues) is 
discarded from the matrices O and Q, and the information of the system carried by the density matrix is reduced. 
However, remaining eigenvectors (if m is large enough) usually describe the system accurately because the truncation 
error e defined as 

e= E ^ (12) 

{discarded} 

is very small (0 < e <C 1). J2{aii} ^ ~ ^^"^ eigenvectors are assumed to be normalized. The matrices O and 
Q enter the linear transformation as projectors mapping two blocks TiWb onto one block T/ through the following 
procedure 

r/^'+^^(ai^[-"'a,+5kiC,'"^"'a;-+5) = E oKa-«'|6a,+4)T/^'Vi6^,+4ki^;a;.+4) 

xW^B(a,+4a,+5k;-+4f^;+5)QKe,''^-+4l^/' "^")- (13) 

Application to the right block T,. is straightforward. As it is seen, we calculate the blocks T/ and Tr separately not 
using the standard mirror-reflection of T; to T,.. This procedure is necessary when dealing with anisotropic and/or 
inhomogeneous systems. 

The calculated new blocks T^^'^^'^ and Tr^^'^^ are used in the next step of the ISM for construction of the new 
superblock 

T[2,+6] = t/^'+i) W^bT^^+i) Wb . (14) 

Within the FSM, e.g., for a sweep to the right only the left blocks are calculated and j'^t^/^^'^) jg taken from the 
previous left sweep 

T[L] = t/^/'-'+'^VbT;(^/'-'-'^W^b. (15) 

The variable k (indexing the steps within a sweep) runs over the values (— fco, — fco + 1, • • ■ , ^ 1, ^o); where 2^° < 
2^/2-2 _ ^YiQ process of sweeping one of the Boltzmann weight is fixed (the upper one in Fig. 1) and the second 

one changes its position within the interval of 2fco lattice sites. The local physical quantities are calculated at the 
lattice sites of the fixed Boltzmann weight and due to the rotational invariance of the problem are valid for all the 
rows of the periodic lattice. 

IV. RESULTS 

It is well known that the DMRG describes better a strip with open boundary conditions than that with the 
periodic boundary conditions [Q because the precision of the largest eigenvalue of the superblock matrix is increasing 
proportionally to m for open boundary conditions while for periodic boundary conditions only as ^yrn. 
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However, if we are not interested in the largest eigenvalue of a finite-strip transfer matrix but in the estimation of 
the free energy of the whole 2D lattice (per spin) , it is more effective to use a strip with periodic boundary conditions 
than that with open boundaries, as demonstrated in Table The results with m = 25 practically exactly reproduce 
the exact values for L ~ 16. The estimation of the free energy for 2D models performed by DMRG can be improved 
by increasing the width of the strip. For a given m the best results are obtained for L — > oo, but in this case, for 
T = 2.1 (i.e. below the critical temperature), the symmetry of the system is spontaneously broken. Exact free energy 

per site /onsTgor "^^^ taken from |l|] . 



FIG. 3. Free energy per site /ism for the Ising model calculated with the standard DMRG method only with the ISM is 
compared with the free energy per site calculated by the modified DMRG algorithm as well as by the exact diagonalization 
method (EDM). N is the order of either the superblock of DMRG or the exact transfer matrix in EDM. 
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The critical temperature and the properties of the infinite 2D system near the critical temperature should be derived 
from finite-size scaling ideas (as the finite- width strip is at criticality for T = only) . 

For calculation of the critical temperature, the phenomenological renormalization approach of Nightingale Jl7[ have 
been used. Here, the scaling properties of the correlation length, found as logarithm of the ratio of two largest 
eigenvalues of the exact or superblock matrix, are exploited. The product of the inverse correlation length Kn and 
the strip width L should not depend on the L at critical temperature Tq{L) 

1. (16) 



(L + 2)Kl+2 

The accuracy of the approximate critical temperature improves with size of the strip in the case of exact diagonal- 
ization. For DMRG calculations this statement is no longer valid, as for very large L the symmetry of the system 
spontaneously breaks, and the phenomenological renormalization is not applicable any more. Thus, for given order 
of approximation m, there exists an optimum value of the strip width L°p*. This can be estimated from the following 
considerations: For exact diagonalization or DMRG calculations with m close to 2^/^~^, the difference of the approx- 
imate critical temperature from the exact critical temperature j'^'^^'^'^^ — 21n^^(l + \/2) [ p^ scales with the width of 
the strip as follows fl^ : 

^c^^^-^c L-''\ (17) 



i.e. the ratio 



rp(cxact) 



R^^E^ = ^L^L. (18) 

The optimum width L°p* should be less than Lq for which the ratio of the derivatives R{Lq) (|iJ) is substantially 
deviated from the originally linear behavior. In our calculations we have considered the DMRG results to be incorrect 
for i? = or oo. In the case of i? = 0, the precise value of L°p' is not too important as the first derivative or change 
of T^{L) is very small. Near R = oo, a sharp drop of the second derivative of T^{L) to zero is required; indeed, 
the change of the distance from the line ^Tpx ^ by more than one order of magnitude takes place within one step of 
strip- width enlargement. 
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FIG. 4. Critical temperatures Tq for Ising model as functions of lattice size L for various sizes of multi-spin variables m 
from DMRG and finite-size scaling. The results for open boundary conditions (OBC) are plotted as triangles while the results 
for the periodic boundary conditions (PBC) are plotted as circles. The exact critical temperature is at the dot-dashed line. The 
OBC plot for m = 30 is indistinguishable from the curve for m = 64 in this figure. 

In Fig. ^ plots of strip- width-dependent critical temperatures Tq{L) for two different boundary conditions and 
various block sizes m are given. The estimations of the exact critical temperature for periodic and open boundary 
conditions were found as the values of Tq(L — 2) if the first or second derivative of T^{L) changed their signs with 

respect to the value in the previous step. The curves for PBC cross the exact value of xhe curve maxima 

for OBC are quite far from it, and by increasing L, Tq approaches the exact value very slowly. 
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FIG. 5. Critical temperatures Tq in Ising model vs. lattice size L and m for the periodic boundary conditions (PBC) only. 
Filled circles represents data which are accepted whereas the open circles are taken as incorrect due to violation of the condition 
(p^. Zoomed Fig. The critical temperature estimations in the inset are given by the rightmost filled circles for respective m. 
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The accuracy of the results for periodic boundary conditions (Fig. ||) is very high aheady at small values of m 
and exceeds by an order the critical temperature estimation for maximum computer-accessible m when using open 
boundary conditions. The critical temperature for not extremely large m = 80 is given to seven digits. As the width 
of the strip can be increased only in discrete steps and the criterion of the optimum strip-width is somewhat vague, 
the accuracy of the critical temperature determination should be taken as large as a single step change of Tq{L). 
These accuracy estimations together with deviations of our results from the exact critical temperature are given in 
Table p. It should be noted that only ISM was performed in calculations of Tq(L) in Figs. ^ and |[ The calculations 
with the FSM has been also done near the L°p' but only slight improvements of critical temperature were obtained. 
In calculation of the thermal critical exponent a similar accuracy was reached; e.g., for m — AA within FSM the 
critical exponent v = 1.0000016 is very close to the exact value v = 1 |p^ . 



FIG. 6. Changes of Tq{L) per one step of strip-width enlargement as well as the deviations of our results from the exact 
critical temperature y^™"'^*' for increasing parameter m. 
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V. CONCLUSION 



The DMRG method for classical spin lattice strips with periodic boundaries was developed and applied to 2D Ising 
model. It was shown that this approach lead to more accurate results for 2D infinite lattice than DMRG with open 
boundary conditions. It was demonstrated that applying finite size scaling to strips treated by DMRG, an optimal 
width of the strip depending on the order of approximation existed, and a prescription how to find Tq(L°p*) was 
given. For the Ising model it was shown by computations that for these L°P'(m) the value of the critical temperature, 
was for a given m closest to the exact one. As our approach does not involve any information about the exact critical 
temperature y^°'^'^'^*) or the universality class of the model, we believe that it is applicable to many different classes 
of spin lattice models. This belief is supported by analogous calculation for anisotropic triangular nearest-neighbor 
Ising model (ATNNI) with two different antiferromagnetic interactions Ji and J2 (model discussed in ||l5[). For this 
model the transfer matrix is non-symmetric and the phase diagram is quite different from that of the standard Ising 
model. For the periodic boundary conditions, the plot of critical temperatures is not monotonously decreasing as in 
the case of the Ising model (Fig. Hj but for large L it turns up. Nevertheless, the accuracy of the critical temperature 
for the exactly solvable case (for external magnetic field H = 0) is similar to the presented ones in this paper. 
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